Charmonium production in relativistic proton-nucleus collisions : What will we learn 

from the negative xp region ? 
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We study the nuclear medium effects on the cc time evolution and charmonium production, in 
a relativistic proton-nucleus collision. In particular, we focus on the fragmentation region of the 
nucleus where the formation length of the charmonium mesons is shorter than the size of the nucleus. 
Little is known on the nuclear effects in this region. We use a quantum-mechanical model which 
includes a realistic potential for the cc system and an imaginary potential to describe the collisions 
of the cc with the nucleons. The imaginary potential introduces a transition amplitude among the 
charmonium states and produces an interference pattern on the charmonium survival probability, 
which is particulartly important for ip' . Our results on the suppression factors are compared with 
data from the NA50 and E866/NuSea Collaborations. Predictions are given for the suppression of 
J/ip, ip' , Xc as a function of the nuclear mass and in the negative xf region, where data will be 
available soon. 



A very exciting era in the study of strong interactions 
has begun with the advent of high energy heavy ion col- 
lision experiments. In particular, this has opened the way 
to the exploration of new states of matter, such as the 
colour glass condensate, corresponding to the nonlinear 
regime of quantum chromodynamics, or the high temper- 
ature phase of nuclear matter, the so-called quark-gluon 
plasma, which is also important for our understanding of 
the early universe. In relativistic heavy ion collisions, 
evidence for the formation of the quark-gluon plasma 
may come from the combined information of different 
signals, such as strangeness, dilepton, photon or charmo- 
nium production. 

Matsui and Satz [1] first suggested that the suppression 
of charmonia could signal that a phase transition has 
occurred. This pioneering idea has triggered a series of 
experiments at the CERN SPS [2,3], whose interpretation 
has been very controversial. It has become clear that, 
to use charmonium as a signal, all possible mechanisms 
which can affect charmonium production need to be well 
understood [4,5]. 

In this context, proton-proton (pp) and proton- nucleus 
(pA) collisions have been intensively studied [2-8]. In 
fact, these processes are used as a reference for the case 
of nucleus-nucleus (AB) collisions where a critical en- 
ergy density may be eventually attained, producing the 
plasma phase. This should affect various observables and 
in particular, it may lead to an extra (anomalous) char- 
monium suppression. There have been indications re- 
cently that the plasma has been produced, but it is still 
too early to draw definite conclusions [9]. These stud- 
ies will be pursued by the running experiments at RHIC, 
and future measurements at LHC. 

It takes a certain time for a cc produced in a collision 
to become a colour singlet state, then it expands to the 
size of a ip meson (the symbol ip stands from now on for 
any charmonium meson, i.e. J/ip, ip' and % c ). Different 



models exist at present for the production mechanisms 
[4,10,11]. 

Sufficiently fast cc can traverse the nucleus before a ip 
meson is fully formed; whereas slow cc can traverse the 
nucleus as fully formed ip meson [12,13]. Experimentally 
one can pinpoint to these different kinematical regions by 
measuring charmonium production at positive and nega- 
tive small Feynman xf and/or making inverse kinematics 
measurements [4,13]. 

One of the aspects which need to be well understood is 
the effect of the nuclear medium on the cc pair evolution. 
In particular, very little is known on the nuclear medium 
effects when the formation time of a fully developed char- 
monium is smaller than the time to traverse the nucleus 
[13] . Soon, measurements in pA and AB collisions will be 
performed by the HERA-B Collaboration [14] at DESY 
and the PHENIX Collaboration at RHIC [15]. Moreover 
the El 60 Collaboration at SLAC will perform similar ex- 
periments where the charmonia are photoproduced [16]. 

In this letter, we focus on this region, i.e. the frag- 
mentation region of the nucleus (negative xf region or 
slow cc) and study the effects of the nuclear medium on 
the cc time evolution and on charmonium suppression 
in relativistic pA collisions [17]. Because of the kine- 
matical region we are interested in, we assume that the 
produced cc becomes quickly a colour singlet state (pre- 
meson) and on its way through the nucleus, this pre- 
meson expands to a ip meson, while experiencing colli- 
sions with the nucleons. To describe this, we use here 
a quantum-mechanical model where, contrary to previ- 
ous works [4,12,13,18], the cc pair is bound by a real- 
istic potential [19] and the premeson (and then meson) 
wave function is expanded on the basis given by the char- 
monium eigenstates. An imaginary potential depending 
on the dipole charmonium-nucleon cross section [12,20] 
is included to describe the collisions with the nucleons. 
This introduces a transition amplitude among the char- 
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monium states. (We assume the premeson wavefunction 
has a fixed angular momentum quantum number, while 
it has different radial components.) 

We will show that the interaction of the premeson and 
later on the meson with the nuclear medium produces an 
interference pattern on the charmonium survival prob- 
ability. We will show how this affects the charmonium 
suppression both when the path in the nucleus is var- 
ied, by changing the mass number A, and as a function 
of xp. These effects are particularly important for ip' . 
We compare our results with experimental data obtained 
at CERN SPS by the Na50 [7] and at Fermilab by the 
E866/NuSea [8] Collaborations. We present predictions 
for the production of J/ip, ip', Xc m a range of small 
negative xp values, for different nuclei [17], for which 
experimental data will be avalaible soon [14,15]. 

Our cc pair, produced in a pA collision, is described 
(in its center of mass frame) by the internal Hamiltonian 

H = ^ + V(r), (1) 
m c 

where m c is the charm quark mass. The potential V(r) 
which binds c and c is a realistic potential which repro- 
duces well the properties of the charmonium family [19]. 
The spin-dependent terms are neglected. Only transi- 
tions between charmonium states with different n quan- 
tum numbers are considered. 

By solving the static Schrodinger equation with Hq 
Eq.(l), one gets the eigenenergies E n ^ and the corre- 
sponding wave functions \n,£) (to simplify the notations, 
we do not write explicitly the dependence on magnetic 
quantum number m). The IS, 2S states correspond in 
our model to J/ip, ip' respectively and IP to \c where 
the different total spin state of % c are not split in en- 
ergy. (We consider that our \ c is produced through a 
two-gluon mechanism). 

The time dependent wave function for the cc in its 
rest frame is expanded on the basis of eigenstates of H a 
Eq.(l): 

OO 

\cc,t)(T) = J2cnAT)e- iE ^ T \n,l). (2) 

(fi, c = 1). In practice, one truncates the sum to a num- 
ber n of eigenstates, for a given i value. As we will dis- 
cuss, we have chosen n large enough so that the results 
we present are not sensitive to this truncation. 

To model the interaction of the cc with the nuclear 
medium we add an imaginary part to Eq.(l) : 

iW = i—(j{f T ,^/s^)p{b,z) (3) 

where v is the speed of the nucleus with respect to the 
cc frame and a is the dipole cross section associated to 
the interaction of the cc with a nucleon N, rp is the 



transverse distance between c and c, ^/~s^n is the en- 
ergy in the center of mass of the ipN system and p 
is the nuclear density evaluated at the position (6, z) 
of the cc center of mass, z being the beam direction. 
For the dipole cross section we use the parametrization 
a(r T ,\fs) = oo(s)(l - e~ r ^/ r "( s )), determined by fitting 
deep inelastic scattering data and which well reproduces 
the charmonium photoproduction data [20]. Concerning 
the nuclear density p, we present results obtained with a 
Woods-Saxon profile, with parameters chosen such as to 
well reproduce the nuclear radii. 

The time evolution of cc wavefunction (2) in the 
nucleus is determined by solving the time dependent 
Schrodinger equation for H = H — iW Eqs.(l) and (3). 
This leads to a system of first-order coupled-channel dif- 
ferential equations for the amplitudes c n j(t) : 

n 

c n , e (t) = -a^CkAty^-^HnJWlkJ) (4) 
k=i 

with a — vp(b, z)/2. The time t — jt is now in the 
laboratory frame. We have neglected in this first calcu- 
lation the higher Fock states that should emerge from the 
Lorentz boost [20]. We see from (4) that the imaginary 
potential (3) introduces a transition amplitude among 
the charmonium eigenstates. 

A difficult choice is that of the initial conditions c ra ^(0) 
for (4), related to the mechanism of hadroproduction of 
charmonia, which is still badly known [4,5,10,11]. We 
use: 

4c-cA0)=c e f e (r T )e-^ 2 ^+^ (5) 

where ct is a normalization constant. We describe the 
conversion of a gluon into a cc through a gaussian multi- 
plied by the function fe(rp) — r\ (or rp) to account for 
the two (one) supplementary gluons necessary to produce 
J/ip,ip' (x c ). The initial wave function depends on one 
parameter only, j3. This parameter has been determined 
by fixing the ratio |c2,o / c i.o| 2 to the experimental ratio 
of ip' over J/4>, produced in pp collisions at 450 GeV, 
i.e. By^&m-**' /Bj/^^aP^ J ^ = (1.60 ± 0.04)% 
[21]. where B^^^^ is the branching ratio to dimuon 
production and a pp ^ the pp reaction cross section. 
This leads to the ratio |c 2 , / c i,o| 2 (0) = 0.22 ± 0.03 and 
(3 = 1.33 ± 0.5 GeV (the same (3 has been used for the 
Xc states). 

The number n of states in (2), to be included in the 
coupled channel equations (4), has to be large enough 
that the results do not depend on the truncation. We 
have checked how the results depend on the inclusion of 
extra states. Up to 4 channels have been used for the S- 
states and up to 2 for the P-states. The truncation always 
affects a little the last state included; whereas the results 
for the lowest energy states are practically unchanged. 
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Let us now come to the results, obtained in infinite 
nuclear matter first, i.e. by taking in (4) a constant nu- 
clear density p = p . In fig.l we show the time evolution 
of the probabilities associated to J/tp and tp' produc- 
tion obtained by solving (4) with the initial conditions 
(5), both neglecting and including the transition ampli- 
tudes between different eigenstates (non-diagonal terms) . 
While in the former case the probabilities decrease expo- 
nentially, in the latter they present an oscillation. The 
interference pattern is more pronounced in the case of tp' 
whereas for J/tp the oscillations stays very close to the 
exponential. This effect directly influences the suppres- 
sion factors in pA collisions as we will see. The dominant 
oscillation frequency in Eq.(4) is uo = (E 2 ,e — E\ t t)/*y. 
We have also seen that the deviation from the uncoupled 
solution becomes stronger the higher 7 is. 

Let us now consider a cc which is produced in a rela- 
tivists pA collision and evolves according to Eqs.(4). In- 
tegrating over all possible paths the ratio |c^,(oo)/c^(0)| 2 , 
which gives the survival probability of the tp in the nu- 
cleus, one gets the pA — > ip reaction cross section : 



a 



.pA—*t/) 



= J dbdzp(b, z)a 



pN — >ip 



c^(t(b,z)) 



£V(0) 



(6) 



where t(b, z) is the time necessary to a cc produced at a 
point (b, z) to traverse the nucleus. If a pA ^ = Aa pN ^ , 
the suppression factor, defined as 



a 



pA — >ip 



AaP N ^ 



(7) 



becomes S\ = 1. 

Before making predictions, let us compare our calcu- 
lations to existing data. For J/tp, we have to take into 
account that the dimuons measured in the detector come 
from the decay of the directly produced J ftp's as well as 
those produced by the decays of i//'s and Xc's. We will 
denote by "J/tp" the total number of produced J/tp, i.e. 



sj'+ = o.62S. J/ * 



0.30S*1 



0.08Sf A [6] 



pA ' pA ' ' pA ' pA 

Let us first discuss charmonium production as a func- 
tion of the nuclear mass. The cc pair travels along dif- 
ferent lenghts, spending different time intervals in the 
nuclei. Therefore, experiments with different A explore 
the time dependence of the c n ^ amplitudes (Fig.l). We 
present results on " J/tp" and tp' (Fig. 2) in comparison 
with the experimental values by the NA50 Collaboration 
[7], obtained for pA collisions with an impinging pro- 
ton energy of 450 GeV. The data are given as branching 
ratio times cross section divided by the nuclear mass, 
which differs from the suppression by a constant. (For 
the comparison, we normalize the data by this constant 
determined from the average of the ratio experimen- 
tal/calculated values.) We can see that our results are 
in good agreement with the experiment both for "J/tp" 
and tp' . In Fig. 3 we give predictions for the suppression 



of Xc as a function of the nuclear mass, for different val- 
ues of 7. The largest suppression is observed when 7 is 
the lowest, corresponding to a longer time spent in the 
nucleus. 

Let us now look at the xp dependence. The sup- 



pression factor is often parametrized as S\(xp, y^Spp) = 
Aai ,(x F ,^,A)-i [4]. i n fig.4 the a values for "J/tp" and 
tp' are compared to experimental data on pA collisions, 
obtained by the E866/NuSea Collaboration [8]. The im- 
pinging proton energy is 800 GeV and the label A is the 
ratio of the atomic numbers of tungsten over beryllium. 
In this region of small positive and negative xp values, 
our calculations are in good agreement with the experi- 
mental values for "J/ tp" , whereas for tp' they overpredict 
the data. In fact, in the region of small positive xp, the 
cc can traverse the nucleus before developing into a tp 
meson. In this kinematical regime, other effects missing 
in our model, like for example gluon shadowing [22], can 
play a role. 

Finally, we present predictions for the a values (Fig. 5) 
and for the suppression factors (Fig. 6) for " J/tp" , tp' and 
Xc in the negative xp region, where experimental data 
from the HERA-B Collaboration [14] at DESY and the 
PHENIX Collaboration at RHIC [15] will soon be avail- 
able. As we can see, both " J/tp" and Xc present a rather 
flat behavior, whereas tp' is very sensitive to the xp val- 
ues. Besides, the results on tp' also show a stronger vari- 
ation with the mass of the nucleus. This can be under- 
stood from the fact that the effect of the coupling between 
the charmonium eigenstates is much stronger for tp' than 
for J/tp (and Xc), as one can see from Fig.l. 

Fig. 5 also shows that the a parametrization of the sup- 
pression (often used with a independent of ^4) is espe- 
cially not good in the negative xp region, since there is 
a strong dependence of the a values on the nuclei. 

In summary, we have studied the production of char- 
monium ( J/tp, tp' , Xc) m relativistic proton-nucleus col- 
lisions. We focussed on the nuclear effects on the cc time 
evolution and on the charmonium production when the 
time for a cc to develop to atp meson is short compared to 
the time necessary to traverse the nucleus. This kinemat- 
ical region may be explored experimentally by looking at 
slow enough cc and therefore charmonium produced in 
the negative xp region. The quantum-mechanical model 
used here includes a realistic potential binding the cc 
pair and an imaginary potential, describing the collision 
of the cc with the nucleons and depending on the dipole 
charmonium-nucleon cross section. The imaginary po- 
tential introduces transition amplitudes among charmo- 
nium eigenstates. We have shown that this produces an 
interference pattern on the charmonium survival proba- 
bility, particularly important for tp' , and affects the char- 
monium suppression factors. The results are compared 
to experimental data from the NA50 and E866/NuSea 
Collaborations. We have made predictions on the sup- 
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pression factors of J/tp, ip' and Xc, for different nuclei 
in the fragmentation region of the nucleus, where experi- 
ments will be performed soon by the HERA-B Collabora- 
tion at DESY and the PHENIX Collaboration at RHIC. 
We expect that the exploration of the negative x p region 
will bring new information on the charmonium-nucleon 
cross section and on the initial conditions which arc not 
well known, especially for \c- A better understanding of 
the mechanisms for charmonium production in proton- 
nucleus collisions in this yet unexplored region will help 
us in the interpretation of the measurements on relativis- 
ts nucleus- nucleus collisions. 
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FIG. 1. Probability of surviving in the nuclear medium 
of a J/tp (upper) and a tp' (lower) as a function of time 
(in the laboratory frame), normalized to the initial proba- 
bility. The transition amplitudes among different charmo- 
nium states, due to the medium, produces an oscillatory be- 
haviour (solid line) whereas the solution is simply an expo- 
nential (dashed line) when they are absent. 
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FIG. 2. Results for the tp production cross section times 
the branching ratio to dimuon production, B^a* , for the 
11 J/tp" (top) and tp' (bottom) produced in pA collisions with 
E p = 450 GeV as a function of the nuclear mass. The ex- 
perimental data obtained by the NA50 Collaboration [7] are 
shown for comparison. (The curves are only given to guide 
the eye) . 
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FIG. 5. Predictions for "J/V>" (top), ?/>' (middle) and X c 
(bottom), on the dependence of the a values as a function of 
xf, obtained with E p = 800 GeV. We present results obtained 
for tungsten (triangles) and beryllium (squares). 
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FIG. 3. Predictions on the suppression factor of the Xc me- 
son, produced in relativistic pA collisions, as a function of the 
mass of the nucleus A. We present results for 7 = 5 (circles), 
16 (squares), 21 (triangles). 
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FIG. 6. Predictions for the suppression factors for u J/ip" 
(circles), ip' (squares) and \c (triangles) as a function of 
the negative xf region, obtained in a pA collision with 
E p = 800 GeV and A being the ratio of tungsten over beryl- 
lium. 
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FIG. 4. Comparison between calculated and experimental 
a values for "J/tp" (top) and tp' (bottom) as a function of 
Feynman xf obtained for pA collisions with E p = 800 GeV. 
Here the a values are relative to the ratio W/Be. The data 
are from the E866/NuSea Collaboration [8]. 




